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Abstract 

This paper suggests a model for a bubble universe using buildable false vacuum bubbles. We 
study the causal structures of collapsing false vacuum bubbles using double-null simulations. 
False vacuum bubbles violate the null energy condition and emit negative energy along the 
outgoing direction through semi-classical effects. If there are a few collapsing false vacuum 
bubbles and they emit negative energy to a certain region, then the region can be approximated 
by a negative energy bath, which means the region is homogeneously filled by negative energy. 
If a false vacuum bubble is generated in the negative energy bath and the tension of the bubble 
effectively becomes negative in the bath, then the bubble can expand and form an inflating 
bubble universe. This scenario uses a set of assumptions different from those in previous studies 
because it does not require tunneling to unbuildable bubbles. 
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1 Introduction 



The information loss problem [T] is related to the issue of whether or not the quantum theory of 
gravity is unitary. It is unclear how we can reconstruct the original information of a black hole after 
the black hole has totally evaporated. If a theory is not unitary, we cannot expect a fundamental 
predictability from the theory. Therefore, many researchers have tried to explain how a black hole 
can conserve information during the evaporation [2] [3] [4] . 

An even worse situation occurs if there is a wormhole or a bubble universe that includes a 
second asymptotic region. If a wormhole metric is given, it will violate unitarity 5 ; and if there is 
a sufficient amount of exotic matter, we can maintain the wormhole for quite long time [6 . However, 
it is unclear how to build such an initial condition. In many cases, the initial condition already 
violates unitarity so that a wormhole is not related to the information loss problem. 

There have been some discussions suggesting that we can derive a bubble universe from a false 
vacuum bubble in a true vacuum background [7 . According to the thin shell approximation, we 
can obtain a bubble universe solution that inflates inside of a Schwarzschild wormhole [8] [5] [10]: 
however, if we assume the null energy condition and global hyperbolicity, whenever we obtain such 
inflating bubble solutions, the initial conditions suffer from the past incompleteness [11] . Therefore, 
unless we violate the null energy condition, wc cannot derive a bubble universe. Of course, if we 
assume arbitrary exotic matter by hand, we can induce a bubble universe, but it is not justified 
unless the matter is allowed by a fundamental theory, e.g., string theory. A better idea to overcome 
the null energy condition was to assume a tunneling, the so-called Farhi-Guth-Guven tunneling 
[12] , The probability of such tunneling was calculated though it was extremely low. However, if it 
is possible in principle, it may imply that the nature allows a violation of unitarity. 

On the other hand, when the geometry of the background is an anti de Sitter space, according 
to the AdS/CFT correspondence [12] (and if we trust the correspondence), the time evolution in 
the anti de Sitter background should be unitary [T3] (this idea may be extended to a de Sitter 
background Then, we can choose between two possibilities: 

1. The generation of a bubble universe does not necessarily imply a violation of unitarity, 

2. Farhi-Guth-Guven tunneling is impossible in such a situation. 

For the first possibility, we can guess that although a Schwarzschild wormhole (or black hole) 
separates two universes and information passes from the first to the second asymptotic region, 
Hawking radiation from the black hole can emit information to the first asymptotic region. This 
should imply a copy of information on a space-like hypersurface. Although such a copy is improb- 
able, if no observer can see the copy of information, we can relax; this is known as the black hole 
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complementarity principle [2J. Although local physics seems to be violated, no one can see the 
violation, and then, it is "a perfect crime." 

However, according to our previous papers, we strongly suspect that black hole complementarity 
is not a correct conjecture |16j|17j. If our opinion is correct, then the reasonable next choice is to 
think that such exotic tunneling should not be allowed in the anti de Sitter background. Apart 
from our opinion, the second possibility is also the consensus between people who are involved in 
holography. Only classically buildable bubbles that do not contain inflating regions are allowed by 
unitary processes |14j . 

In this controversial context, we will argue that a bubble universe solution is possible without 
assuming unbuildable bubbles. We will only assume buildable bubbles, where the term 'buildable' 
means that the bubble does not contain an inflating region. However, we can violate the null energy 
condition through Hawking radiation of the de Sitter space |18) ; a false vacuum bubble will emit 
Hawking radiation and will violate the null energy condition [19 20 . Then an outside observer of 
the bubble can obtain a small region that can be approximated by a negative energy bath, which 
means the region is homogeneously filled by negative energy. If we further assume that there is 
a second tunneling of a false vacuum bubble in the negative energy bath, then the region allows 
expanding and inflating bubble solutions via a violation of the null energy condition, even though 
the bubble is classically buildable. 

In this paper, we will confirm this idea using numerical simulations. We followed the double- 
null formalism which was developed by previous researchers [T7] [2T] [32] . In particular, to include 
Hawking radiation of a de Sitter space, we assumed S'-wave approximation using 1 + 1 dimensional 
results ;19 [23 . This approximation gave good results for black holes [T7] | 21 j . 

Section [2] summarizes previous results on dynamics of false vacuum bubbles in the thin shell 
approximation. In Section [31 we develop a model that is beyond the thin shell approximation for 
a false vacuum bubble and includes semi-classical effects. Then we discuss the generation of a 
negative energy bath using false vacuum bubbles. In Section |4j we discuss a generation of a bubble 
universe using the negative energy bath and in Section [SJ we discuss on our assumptions for bubble 
universes. Section [5] summarizes and contrasts the ideas already known and our new contribu- 
tions. In Appendix A, we check the consistency and convergence of our numerical calculations; in 
Appendix B, we compare our numerical setup to analytic expectations. 
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Figure 1: Solutions of the thin shell approximation (symmetric cases). Since $ is always positive 
in the r — > limit, cLSd is disallowed; /?+ is always positive in the r — > limit, Schc is disallowed; 
is always negative in the r — > oo limit, Schs is disallowed. 
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Figure 2: Solutions of the thin shell approximation (asymmetric cases). 




2 Dynamics of false vacuum bubbles in the thin shell ap- 
proximation 

If we can apply the thin shell approximation [24] with spherical symmetry, the inside can be de- 
scribed by a static dc Sitter space and the outside by a static Schwarzschild space: for the inside, 

ds 2 _ = -f-(R)dT 2 + —l-—dR 2 + R 2 dn 2 , (1) 



and for the outside, 



f-(R) 



ds\ = -f+(R)dT 2 + — — j—rdR 2 + R 2 dn 2 , (2) 



f+(R) 

where f-(R) = 1 — R 2 jl 2 and f+(R) = 1 — 2M/R. The remaining part is on the dynamics of the 
shell at R = r(r) and the dynamics is determined by 



e-y/r 2 + f--e +y /r*+f + = Aura, (3) 
r 2 + V cff {r) = 0, (4) 



where 



Ven(r) = U , (5) 

t is the time parameter of the shell, a is the tension of the shell and e± are +1 if the outward 
normal to the shell is pointing towards an increasing r and —1 if towards a decreasing r [7] [8] [5] [13] . 
In general, the effective potential V c s is a convex function for a time-like shell |8] [5] [H] ; therefore, 
it allows a collapsing solution or an expanding solution. There are basically five possibilities: (a) 
from expanding to collapsing, (b) from collapsing to expanding, (c) from collapsing to collapsing, 
(d) from expanding to expanding, and (e) a static solution in an unstable equilibrium. Solutions 
(a) and (b) are symmetric solutions, whereas (c) and (d) are asymmetric solutions. For simplicity, 
we omit the unstable equilibrium case (e). 

To maintain the information of the sign of each of the roots e± , we need to compare the extrinsic 
curvatures for the outside and the inside of the shell [5] [5] Q3] . The extrinsic curvatures are defined 
as follows: 

2„2J1 



P- = 5 = ±V r + J-, (6 

8nra 



and 



/- - /+ - 167r 2 cr 2 r 2 



8 + = J J \ = ±V^TU- (7) 

8nra 
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Now, to satisfy the Einstein equations, the relation 

/3_ - 13+ = Anar (8) 

should hold Q3]. 

When the null energy condition is satisfied, we can classify all causal structures and some 
solutions show formation of bubble universes (Figures Q] and 0) [5] [5] [H] . 

First let us classify symmetric solutions. The left diagram of Figure [T] is for the de Sitter space, 
and the right diagram is for the Schwarzschild space. For a collapsing case, dSA or dSu are possible; 
and Sche, Schc or ScItd are possible. For an expanding case, dSe or dSc are also possible; and 
ScIia or ScIie are possible. However, according to the behavior of the extrinsic curvatures in the 
r — > or r — > oo limit, we can remove the solutions of dSn, Schc and ScIie- Therefore, there are 
four possible solutions: dSA — Sche, dSA — ScIid, dSe — ScIia and dSc — ScIia- Case dSA — Scire 
is a collapsing bubble solution, where the collapsing shell is inside of a Schwarzschild wormhole. 
Case dSA — Scho is a collapsing bubble solution, where the collapsing shell induces a Schwarzschild 
black hole. Case dSe — ScIia is an expanding bubble solution, where the shell expands inside of 
a Schwarzschild wormhole and the shell becomes greater than the horizon size of the inside de 
Sitter space. Case dSc — ScIia is an expanding bubble solution, where the shell expands inside of 
a Schwarzschild wormhole and the shell expands outside of the cosmological horizon for the r = 
observer. 

Second, let us classify asymmetric solutions (Figure^]). The most interesting case is the creation 
of a bubble universe. In this case, we need to consider the expanding to expanding solution. Here, 
dSE, dSF and Schp are allowed, thus giving us the case of dSE — ScIie and dSp — Schp as allowed 
transition solutions. We can interpret these as expanding solutions that begin from a singularity. 

Note that the initial condition should be singular for such bubble universes, if we assume general 
relativity, global hyperbolicity and the null energy condition [11] . Hence, if a false vacuum bubble 
includes inflation in its interior, its initial condition is unbuildable in the general relativistic sense. 
Therefore, to generate a bubble universe, some previous researchers considered tunneling from a 
classically buildable bubble to a classically unbuildable bubble (Figure [3]) and the tunneling rate 
could be calculated [T2] . 
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3 Generation of a negative energy bath using false vacuum 
bubbles 



In this section, we study the dynamics of collapsing false vacuum bubbles. In the thin shell ap- 
proximation, the dynamics was already studied in the previous section. However, in our analysis, 
we try two different points: (1) we go beyond the thin shell approximation and (2) we include the 
semi-classical effects on the false vacuum bubble. 

3.1 Setup 

We describe a Lagrangian with a scalar held $ and potential V($) [2"5][2"6]: 

C = -^, a ^, b g ab -2V{^). (9) 
From this Lagrangian we can derive the equations of motion for the scalar held: 

<S>:ab9 ab -V'{<P) = 0. (10) 

In addition, the energy-momentum tensors become 

T ab = $;a$;b - \g a b (® ;c$;d9° d + 2V(*)). (11) 

Now, we describe our numerical setup. We start from the double-null coordinates (our convention 

is [u,v,0,<p]), 

ds 2 = -a 2 (u, v)dudv + r 2 (u, v)dfl 2 , (12) 

assuming spherical symmetry. Here u is the ingoing null direction and v is the outgoing null 
direction. 

We define the main functions as follows (we follow the numerical approach of previous authors 
[27] [17] [21] [22].): the metric function a, the area function r, and the massless scalar field S = VAttQ. 
We also use some conventions: d = a. v /a, h = a. u /a, f = r tU , g = r „, W = 5 U , Z = S lV . 

From this setup, the following components can be calculated: 

Gun = --(/,« -2fh), (13) 
r 

Guv = ^ (4r/^ + a 2 + 4/ 5 ) , (14) 

G vv = --(g. v -2gd), (15) 
r 

Ggg - -4*(d t „ + &), (16) 
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47T 



W 2 , 



T 
T 



±z 2 , 

4tt ' 



2irot- 



■WZ~r 2 V{S), 



where 



v(s)=vm^ s/v ^. 

From the equation of the scalar field, we get the following equation: 

rZ tU + fZ + gW + Tra 2 rV' (S) = 0. 



(17) 
(18) 
(19) 
(20) 

(21) 

(22) 



Note that, V (S) = dV(S)/dS. 

We also consider renormalized energy-momentum tensors to include semiclassical effects. Spher- 
ical symmetry makes it reasonable to use the 1 + 1-dimensional results [53] divided by Airr 2 [T7] [2T] : 

P 



Anr 2 



(h u - h 2 ) 



P 

Anr 2 



(d v - d 2 ) 



(23) 
(24) 
(25) 



with P = Nlp^/12-jr, where N is the number of massless scalar fields and lp\ is the Planck length. 



We use the semi-classical Einstein equation, 

G^ v = 8"7T (j 1 ^ + (f^ u )^j . 
Finally, we summarize our simulation equations: 
1. Einstein equations: 



(26) 



1 



wz 



p 



2dg-rZ 2 (d v -<T), 

r 



g,u = f.v = -— - t~ + 2na 2 rV(S) ~ —d, u , 
r 4r r 



f, u = 2fh-rW 2 



P 



2. Scalar field equations: 



fZ gW 



-(h u ~h 2 ). 



ira 2 V (S). 



(27) 

(28) 
(29) 
(30) 

(31) 
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3.2 Initial conditions 

First, we use the following potential: 

S < 0, 

i-(S')= ^ A[i_cos(^)] 0<S<oj, (32) 
A w < S. 

We prepare a false vacuum bubble along the initial outgoing surface. We need the initial 
conditions for each function on the initial u = Ui = and v = u; = surfaces. There are two kinds 
of information: geometry (a, r, g, /, ft, d) and matter (S, W, Z). 

On the geometry side, we have gauge freedom to choose the initial r function; although all 
constant u and v lines are null, there remains freedom to choose the distances between the null lines. 
We choose r(u, v{) = ur U Q + tq and r{u\, v) — vr V Q + rg. Here, we fix tq — 10. Then, g(u- u v) = r v o 
and f(u,Vi) = r U Q are naturally obtained. For convenience, we choose r U Q = —1/2 and r V Q = 1/2. 
In general, the metric function a(u, Vi) has a smaller value than that of a pure de Sitter space if we 
study the inside of a false vacuum bubble (see the discussion in Section 4.2.3 in [22])- I n terms of the 
mass function m(u, v) = (r/2)(l + 4r u r v /a 2 - 8ttW 2 /3) [25], a(ui,«i) = (1 - SirVrl/Z)- 1 / 2 > 1 is 
the choice for a pure de Sitter space, and hence it is reasonable to choose a(0, 0) = 1 to study the 
inside of a false vacuum bubble. 

On the matter side, first we fix S(u,v{) = So- Then S(ui,v) will be defined by 

V < Vsholl, 

v) = { (S + S-)G(v) - S- v shen < v < vshdi + A^, ( 33 ) 

w sh cii + Aw < v, 




where 

G(v) = 1 - sin 2 



ir(v- Ushoii) 



(34) 



2Av 

and where SL is a constant. We fix Sq = u), v s hcii = 10 and Aw = 0.5 for convenience. Then, we can 
calculate S(u,Vi), W(u,Vi), S(ui,v) and Z(ui,v) for the initial surfaces. Then, as we fix S(u{, v), 
from the Einstein equations, we can obtain u(ui,v) from 2gd — rZ 2 (since r >vv = at the initial 
surface and asymptotically we can assume that P = 0). Finally, the other functions can be evolved 
using equations on a >uv , r iUU or r <vv and S iUV . 

The remaining free parameters are P (the strength of Hawking radiation), So and S- (the 
amplitudes of the scalar field) and A (the vacuum energy of the inside false vacuum region). 

We can choose the evolution equation among Equations (12"51) , (|2l?|) or (15U1) to obtain the function 
r. In this paper, we used Equation (|29l) . Then the remaining equations become constraint equations. 
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We obtained r by Equation (|28p and ([50)1 and compared it with the original result to check the 
consistency of the simulations (Appendix [2]) • In this paper, the numerical code was described by 
the 2nd order Runge-Kutta method [29 . Its convergence was tested in Appendix El 

3.3 Results 

We report the simulation results in Figures HI [51 [5] and [7J The parameters are P = 0.1, So = 0.01, 
S- = 0.1 and A = 0.0008, 0.00077, 0.00075 and 0.0006. 

In Figure [31 we can see two horizons: r „ = and r_ u = 0. The former is the apparent horizon 
due to the collapse of the shell and there is a space-like singularity inside of the horizon. The latter 
is the horizon of the false vacuum bubble. As the shell collapses, the r M > region disappears. 

If the size of the shell is smaller than the horizon size of the false vacuum bubble, we will not 
see the r tU = horizon (Figure [7]). Between the two limits (Figures @] and [7J , one would expect to 
see a smooth transition: Figures [5] and [5] show the proper transition. If there is no violation of the 
null energy condition, r. u should decrease for all ingoing observers since 

at the r >u = horizon. Hence, if G uu > 0, then r iU cannot increase. Then it is difficult to combine 
the two limits in Figures H and [JJ However, in our setup, we introduced a violation of the null 
energy condition and the violation smoothly connects the two limits (Figure [8]). Near the initial 
surface, we can see a piece of r u = horizons in Figure [5] This is also a remainder during the 
transition, but it is negligible in the whole causal structure. 
Figures l9l fTUl and [TT1 are interpretations of the results. 

For Figure [HI the shell is already larger than the horizon size of the false vacuum bubble and 
it will eventually form the future infinity. As the shell collapses, it will form a black hole and the 
outside observer will see just a true vacuum region. Then to join the inside and outside, an r > 
region appears and disappears. There should be a region where the singularity and the future 
infinity meet; in that region, some complicated things happen, so there should be a Cauchy horizon 
to connect the singularity and the future infinity. This diagram is intuitively similar to dSc — ScIia- 

For Figure[Tl] the shell is sufficiently small so that we cannot see the r „ = horizons. Therefore, 
this is a buildable solution. This diagram is intuitively similar to dSA — ScIld- 

Between the two diagrams (Figures and [IT]), there can be Figure HOI due to the violation of the 
null energy condition. The horizon dynamics occurs on the shell and hence it cannot be described by 
the thin shell approximation. Therefore, this is a new diagram beyond the thin shell approximation 
and demonstrates the semi-classical effects. 
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Figure 8: As A decreases or the size of the shell becomes smaller, if there is no violation of the null 
energy condition, two r„ = horizons will not meet (left). However, if there is a violation of the 
null energy condition, two horizons can be combined (right). 

If we see the energy- momentum tensors T uu and T vv , we can easily notice that T uu < and 
T vv > inside of the shell. This implies that the false vacuum bubble emits negative energy 
along the outgoing direction and emits positive energy along the ingoing direction; the situation is 
opposite for a black hole that emits negative energy along the ingoing direction and emits positive 
energy along the outgoing direction. As the shell collapses, the outgoing negative energy escapes 
to the outside of the shell. 

One interesting observation is that the collapsing shell initially emits positive energy along the 
outgoing direction. In the thin shell approximation, there is no such effect and this effect is not 
an essential one. Even with a thick shell, it is possible to control such effects (Type 1 in |22|). 
In our setup, we may give too much energy to the shell and hence some part of the energy would 
be emitted to the outgoing direction. As we choose sufficiently small energy on the shell, the 
positive energy along the outgoing direction will be controlled. In Figure Q21 we tested some cases: 
P = 0.1, A = 0.0006 and varying S = £_ = 0.001, 0.0005, and 0.00001. We also tested the P = l, 
So = S- = 0.0005 and A = 0.0006 case. Since the energy on the shell T vv ~ ((5 + S-)/Av) 2 
decreases, the black hole horizon r iV — disappeared in the integrated domain. As the energy 
on the shell decreased, the outgoing positive energy decreased. In addition, as we assume strong 
Hawking radiation, the outgoing positive energy decreases. Therefore, although there is a positive 
energy flux along the outgoing direction, it can be controlled and it is not the essential property of 
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Figure 12: T uu for various cases. Since the energy of the shell is sufficiently small, we cannot see 
r „ = horizons in these figures. 

a collapsing bubble. 

In conclusion, it is clear that a collapsing false vacuum bubble emits negative energy along the 
outgoing direction and the negative energy can be observed outside of the shell. This is due to the 
semi-classical effects of the de Sitter space. 



3.4 Generation of a negative energy bath 

To include the semi-classical effect in false vacuum bubbles, we used the S'-wave approximation of 
the renormalized energy-momentum tensor of 1 + 1 dimensions (metric is ds 2 = —a 2 (u,v)dudv) 

MM- 



y^vv) 



p_ ( 

P . 

x [aa . 

or 

P ( 

or 



■ 2a „ 2 ) , 



(36) 
(37) 
(38) 



where P is a constant that is proportional to the number of independent modes of Hawking radiation. 
For a de Sitter space, (T uu ) = (T vv ) ~ —PA |19j . If we assume the S- wave approximation, we will 
obtain (T uu ) — (T vv ) ~ —PA/r 2 in spherically symmetric 3 + 1 dimensions. Then, we can see the 
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violation of the null energy condition from the false vacuum bubbles. 

This possibility was confirmed by our simulations. Figures El [5] El an d HH as well as their 
interpretations, Figures [9j [10] and [TTJ show that negative energy is emitted along the outgoing 
null direction or a violation of the null energy condition for T uu components. As A increases, the 
negative energy increases in Figures 0] [SJ [B] and [7J In Figure [T21 as P increases, the amount of 
negative energy increases. 

As A and P increases, we obtain large violations of the null energy condition. However, A and 
P values are limited. First, the size of the shell should be smaller than I oc 1/y/A. Therefore, A has 
an upper bound: 

A<-^- (39) 

' shell 

Second, a large P gives the cutoff ~ \fP in Equation (|2"7|) . which is a boundary of the semi-classical 
description [50] : 

P < dell- (40) 

Therefore, the amount of negative energy should be bounded: 

\( T )\ZP^<^- (41) 

If we observe a bubble at a sufficiently large distance r from the bubble, then the emitted negative 
energy will be on the order of ~ 1/r 2 . One note of caution is that even though we apply Dvali's 
cutoff 30 , in fact, we do not have a singularity in the center of the false vacuum region. Therefore, 
even though we use a sufficiently large P, if the center does not have a singularity and the final 
observing point is larger than VP~, then the result may still be semi-classical. Thus, it might be 
possible to obtain a lot of energy from bubbles through a large P. 

There is an outgoing positive energy 'bounce' which may annihilate such negative energy flow. In 
our setup, we assumed an ingoing null shell with arbitrary energy on the shell. However, there should 
be a more proper description using a thin shell with a proper tension: dSA — Scho • Figures [7J and 
1111 correspond to the collapsing part of such a solution. For such a proper limit (and in Figure 1121) , 
the 'bounce' will not annihilate the negative energy. 

Now, let us imagine the following situation. In a background, many buildable (i.e., does not 
inflate) false vacuum bubbles are generated and annihilated by collapsing. Then such bubbles will 
emit negative energy along the outgoing direction. At a distance r, the emitted energy will be 
on the order of 1/r 2 . Then, accidently, there can be a certain small region where negative energy 
of various directions are concentrated (Figure IT5|) . In that region, negative energy flows along all 
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Figure 13: In the upper figure, buildable false vacuum bubbles (blue dashed circles) generate 
negative energy along their outgoing directions. Then the energy can be concentrated in a certain 
region (black dashed circle). If we look at the region locally, it looks like a negative energy bath. 
In the lower figure, we consider the dynamics of a false vacuum bubble in the negative energy bath. 
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directions and appears homogeneous: it seems that there is a bath of negative energy photons. In 
this region, we can use the energy-momentum tensor 

Tfm = diag[-p, -p, -p, -p] . (42) 

If the amount of p and p are not too large so that we can ignore the back-reaction to the background, 
we can still assume that the background is nearly flat. We call this background a negative energy 
bath. 

If a false vacuum bubble is generated and the tension of the bubble is sufficiently small so that 
the negative energy background can make the tension effectively negative, then the subsequent 
evolution will be totally different from the original result of pure Einstein gravity. 



4 Generation of a bubble universe using a negative energy 
bath 

4.1 The initial setup of a negative energy bath and a bubble 

Let us assume the energy- momentum tensors = T^J + T$ : 

-j, ! —Ag^ u in false vacuum regions, 

-^jUf \ (^) 

I in true vacuum regions, 

with the negative energy bath 

Tffi = diag[— p, —p, —p, —p\. (44) 

Note that, approximately, the negative energy bath contributes to the tension of the shell by 
T^v e ~ — pe, where e <C l,M (or, M = 0) is the thickness of the shell. If \pe\ > a, then the effective 
tension a = a — pe becomes negative. 

As long as the shell is in the negative energy bath, such an approximation holds effectively. 
However, if the shell expands further than the region of the negative energy bath, the tension 
will be positive again. Therefore, to obtain inflation before the shell reaches the boundary of the 
negative energy bath, we require that the horizon size of the inside de Sitter space is smaller than 
the size of the negative energy bath. Of course, it can be chosen as a free parameter [3T] . 

4.2 Dynamics of thin shell bubbles and generation of a bubble universe 

Here, we consider bubbles with M = limit, i.e., bubbles in a nearly flat background. If A ~ 
1/Z 2 ^ p,p, then the inside and the outside will follow the de Sitter space and the flat space. Now 
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Figure 14: An example of V e g, for I = 2 and a — 0.01. In general, the effective potentials for M = 
are qualitatively the same. 




Figure 15: If the tension is positive, a false vacuum bubble in a flat background is unstable and 
hence buildable. 
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flat 



de Sitter 



Figure 16: Causal structure of a stable expanding bubble. dSe — FlatA and dSc — FlatA are both 
possible. 




25 



the equation of motion becomes [5] [5] [T3] : 



e-y/P + f- -e +y /r 2 + f+ = Airra, (45) 
r 2 + V cff (r) = 0, (46) 



where 



^ff(r-) = /+ (47) 

and /_ = 1 — r 2 /I 2 , f + = 1, r is the radius of the bubble and e± determines the direction of the 
shell. 

To maintain the information of signs of roots e±, we use the extrinsic curvatures again: 

2=2^2 



/_-/+ + le^v 



f>- = - ^T"' - ±V^Tj:, (48) 

and 

/_-/+- 16ttW 



P+ = ^ = ±V^+U- (49) 

Whether a is positive or negative, the effective potential V c g (r) will have the same form as long 
as \a\ and I are the same. In this limit, all bubbles should expand (Figure [T4")) . If a is positive, /3+ 
is always negative and hence it should touch the left boundary of the Minkowski space; or if the 
bubble is in the right part of the Penrose diagram, the bubble should be unstable (Figure ITS"]) [22] . 
Therefore, such a bubble is buildable and can be generated initially. If a is negative, then j3 + is 
always positive and hence it can touch the right boundary (Figure [TBI . Therefore, the bubble can 
inflate and can form a bubble universe (Figure IT7l) . In this case, even though the bubble expands 
and inflates, the initial state of the bubble is not unbuildable. Rather, in this case, a buildable 
bubble became an inflating bubble through a violation of the null energy condition. 

If the amount of negative energy is sufficiently smaller than the vacuum energy of the inside 
region, it will not affect the inside of the shell; if it is sufficiently small, it will not affect the 
outside either; and if the energy-momentum tensor component of the shell is sufficiently smaller 
than the density of the negative energy, the shell can have a negative tension successively. In fact, 
it was confirmed that such a situation (negative energy is concentrated on the shell and the tension 
becomes negative) allows inflation by Type 3 in |22j . 



5 Remarks of caution on bubble universes 

If there are r u — horizons, it is not so strange that there is radiation along the ingoing or outgoing 
directions. However, what will happen if a bubble is slightly smaller than size of the horizon so that 
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one cannot see the horizon? Our calculations confirmed that false vacuum bubbles emit a negative 
energy flux along the outgoing direction not only as seen in Figures and [TU] but also in Figure [TT] 
and dSA — ScIid- In the latter cases, the bubbles themselves do not contain inflation and hence 
they are in principle buildable. 

If this is true, we can obtain a negative energy bath using buildable bubbles (Figure IT31 . In 
the negative energy bath, a buildable bubble can expand and inflate through a violation of the null 
energy condition (Figures 1161 and IT71) . If such bubble universes are in principle possible, then it will 
have theoretical importance for unitarity and holography. 

However, such possibilities of creating bubble universes are very rare and accidental. We cau- 
tiously remark what are really necessary conditions to obtain such negative energy baths and bubble 
universes. 

1. Although our false vacuum bubbles are buildable in the sense that there is no inflating region 
inside of them, it is still unclear whether such buildable false vacuum bubbles can be obtained 
via quantum tunneling. No one can give a definite answer, but if bounce solutions in the 
Euclidean signatures are possible, it will be more probable. As we know, there is no known 
false vacuum bubble bounce solution if the bubble is smaller than the size of the horizon of 
the background de Sitter space [32] [9] . 

However, in scalar-tensor gravity, there are some bounce solutions of false vacuum bubbles 
[3"3"] . and in some cases, the scalar-tensor gravity will be embedded in string theory [34] . 
Therefore, although it is rare, string theory may allow such situations. 

2. In many cases, we observed not only outgoing negative energy fluxes but also outgoing positive 
energy fluxes. Of course, in the thin shell limit, we will control such excessive energy to be 
sufficiently small. However, in realistic and natural situations, the positive energy flux may 
annihilate negative energy baths. 

3. Even though such a negative energy bath is formed, the amount of negative energy will in 
general be very small. Then, to see the formation of a bubble universe, the tension of the 
bubble should be extremely small while the inside vacuum energy remains at a constant value. 
It is not impossible in principle, as long as we tune the potential so that the difference of fields 
between the two vacuums to be sufficiently small [31] . However, it needs a fine-tuning on the 
potential. If a potential has fundamental constraints, then such small tensions will not be 
obtained and hence the negative energy bath will be disabled by large tensions. 

4. In addition, although such tension is sufficiently small, as the shell expands, if the shell 
becomes larger than the region of the negative energy bath, the tension will be positive again. 
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Therefore, the effective cosmological constant of the inside de Sitter space should be sufficiently 
large so that inflation begins before the shell becomes larger than the negative energy bath. 
However, if we assume an arbitrarily fine-tuned potential (e.g., landscape [35]), such small 
tension is possible in principle. 

Therefore, the scenario of the formation of a bubble universe using a negative energy bath 
requires favorable chances and fine-tunings of the potential. However, if it is in principle possible, 
for example in string theory, then it will still have theoretical importance. 

6 Discussion 

We studied the dynamics of collapsing false vacuum bubbles with semi-classical effects beyond 
the thin shell approximation using double-null simulations. Using this setup, we investigated the 
possibility of the generation of a negative energy bath using buildable bubbles. We also studied the 
creation of a bubble universe in a negative energy bath. 

There is a novel comment on the causal structure of a collapsing false vacuum bubble. An 
ingoing observer cannot see r iUU > at the r „ = horizon if there is no violation of the null energy 
condition. If a false vacuum bubble which is slightly larger than the size of the horizon of the 
internal de Sitter space collapses, the r„ > region should disappear, and it seems that r uu < 
at the r u = horizon. As we violate the null energy condition by renormalized energy-momentum 
tensors, we can smoothly connect a de Sitter space and a flat space. Therefore, the semi-classical 
effects regularize and smoothly connect the dynamics of false vacuum bubbles. 

Before concluding this paper, we summarize the already known ideas and our new conclusions 
and suggestions. The following ideas are already known by previous researchers: 

(a) If tunneling from a buildable bubble to a unbuildable bubble is possible, it will violate unitarity 
[12] ; but it is controversial whether or not such tunneling is allowed 14]. 

(b) False vacuum bubbles violate the null energy condition and emit negative energy to an out- 
going null direction [15] [ 15 ] [20] . 

(c) If a violation of the null energy condition is allowed so that the tension of a false vacuum 
bubble can be negative, then an expanding and inflating bubble is buildable and possible 

iguana- 

Now the authors contribute the follow points in this paper: 
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• We numerically confirm the emission of the negative energy of false vacuum bubbles, even 
though the false vacuum bubble is buildable. As we tune initial conditions, the amount of 
negative energy is sufficiently controllable. 

• The authors connect ideas (b) and (c) above to make a bubble universe, without assuming 
unbuildable bubbles: if negative energy (emitted by buildable false vacuum bubbles) is con- 
centrated in a region (idea (b)), then it will be possible to see an expanding and inflating 
bubble in the region (idea (c)). 

First of all, we reduced one condition for a bubble universe to other simple conditions. Previously, 
researchers thought that to obtain a bubble universe from buildable bubbles, we needed Farhi-Guth- 
Guven tunneling. Now we argue that Farhi-Guth-Guven tunneling is not a necessary condition and 
even though we have buildable bubbles, they violate the null energy condition and it will be possible 
to form a bubble universe. The other necessary conditions of our scenario are some fine-tunings of 
the potential. Although it is unnatural, if string theory allows such fine-tunings, this can be a good 
starting point for interesting discussions on information and unitarity. 
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A Consistency and convergence checks 

In this appendix, we discuss convergence and consistency tests for our simulations. 

We mainly obtain the function r by integrating Equation f|2Q|) . However, we can also obtain 
the function r by integrating Equation (|28p or (j3"0")) . To check the consistency, we used another 
scheme: obtain g from Equation (|28p and obtain / from Equation (|3"0|) . We call this result A 2 \ We 
compared r and r*- 2 -* to check consistency for u = 6, 12 and 18 slices. We used parameters P = 0.1, 
So = 0.01, S- = 0.1 and A = 0.00075. Figure [18] shows that they coincide well and the error is less 
than a few percents except near the singularity or deep inside of the black hole. 



29 



10 20 30 40 50 60 

v 



Figure 18: Comparison between two different integration schemes r and A 2 \ where the former is 
integrated by Equation (|2T))) while the latter is integrated by Equation (|28p and Equation ([3U1) . Note 
that two schemes are equivalent except near the singularity or deep inside of the black hole. The 
errors are less than 1% for almost all integrated domains. 

Also, we check the convergence by comparing 1 x 1, 2 x 2 finer and 4x4 finer simulations for 
u = 6, 12 and 18 slices. We used parameters P = 0.1, S = 0.01, S_ = 0.1 and A = 0.00075. 
Figure [19] shows that they converge to the second order and the error is less than 1 % for almost 
all integrated domains. 

B Analytic test for numerical setup 

To trust whether our numerical assumption for the renormalized energy-momentum tensor [19 [23 
gives correct physics, a comparison between the analytic and numerical results will be a good test. 
We know the exact 2-dimensional results (up to 1-loop order) and hence Equations ((23"]) . (j2~4|) . and 
(f2"5|) are the most intuitive and reasonable approximation to the 4-dimensional cases. For black hole 
cases, there are studies confirm that this approximation gives correct thermodynamics |21) . 

Here, we applied this method to a false vacuum bubble. However, there are some obstacles 
to comparing the analytic results and the numerical results in the false vacuum background. The 
known semi-classical results are temperature, entropy, the renormalized energy-momentum tensor 
components, etc., for the static limit or for the de Sitter space. The relevant quantity is to compare 
the renormalized energy-momentum tensor components of the static limit to our numerical results. 
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Figure 19: Convergence test. We can see the second order convergence. Errors are less than 1% for 
almost all integrated domains. 

Note that, in our setup, we prepare a classical background without Hawking radiation and then 
we turn on a field combination or Hawking radiation on the background. (Of course, the Einstein 
equations and field equations are satisfied.) Therefore, energy-momentum tensor components can- 
not be static and they are essentially time dependent. Moreover, the energy-momentum tensors will 
give back-reactions to the metric components and these back-reactions will make the comparison 
difficult. 

In this appendix, we compare the energy-momentum tensor components (T^ u ) as we vary A. 
We have to determine a space-time position where we make comparison. This is not so well-defined, 
but one reasonable position is to see the maximum value of (T^ u ) along an almost asymptotic out- 
going null surface (u ~ 0.0067), since that is relatively free from the back- reactions of renormalized 
energy-momentum tensors. 

Figure [501 is a plot for the maximum of | (T^ u ) | along an asymptotic surface (u ~ 0.0067) 
as varying A = 0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.00075,0.00077,0.0008,0.0009. The 
gradient is 0.86789 ±0.01333 and approximately 1. Note that the smallest value of (T^ u ) is observed 
on an almost same radius. We analytically expect that the renormalized energy-momentum tensors 
in a de Sitter background should be [13] 

(T*) « (5°) 

where ro is a certain constant radius. Therefore, we expect that the proportionality will be approx- 
imately 1. This correspondence is not so trivial since our numerical setup does not contain a direct 
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Figure 20: A log-log plot of \(f^ u ) \ and A. The gradient is 0.86789 ± 0.01333 and approximately 1. 



relation between A and (T^ u ). Therefore, we can conclude that our numerical calculations exhibit 
this behavior consistently, even though there are some ambiguities. 
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